Measuring equilibrium properties in aging systems 
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We corroborate the idea of a close connection between replica symmetry breaking and aging in the 
linear response function for a large class of finite-dimensional systems with short-range interactions. 
In these system, characterized by a continuity condition with respect to weak random perturbations 
of the Hamiltonian, the "fluctuation dissipation ratio" in off-equilibrium dynamics should be equal 
to the static cumulative distribution function of the overlaps. This allows for an experimental 
measurement of the equilibrium order parameter function. 
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The glassy state of matter can appear in systems 
with quenched disorder (like spin-glasses), or in non- 
disordered systems. Ergodicity breaking takes a special 
form in these systems. A rather generic situation is the 
existence of many solid, "glass" , phases, which are very 
different from one another, and unrelated among them- 
selves by symmetry transformations. Hence the Gibbs 
equilibrium measure decomposes into a mixture of many 
pure states. This phenomenon was first studied in detail 
in the mean field theory of spin glasses, where it received 
the name of replica-symmetry breaking JlJ . But it can be 
defined in a straightforward way and easily extended to 
other systems, by considering an order parameter func- 
tion, the overlap distribution function. This function 
measures the probability that two configurations of the 
system, picked up independently with the Gibbs mea- 
sure, lie at a given distance from each other ||. Replica- 
symmetry breaking is made manifest when this function 
is nontrivial. 

The existence of non trivial overlap distributions, 
first found in mean field systems, has been shown un- 
ambiguously, through numerical simulations, in finite- 
dimensional spin-glass systems with short-range inter- 
actions This order parameter function is a very 
important tool for the mathematical description of the 
Gibbs state. Unfortunately it seems impossible to ac- 
cess it experimentally for two reasons: (I) Large glassy 
systems never reach equilibrium at low temperatures; 
(2) The measurement of the distance between configura- 



tions requires a detailed observation at the microscopic — 
atomic — level, which is impossible. (In simulations, the 
second objection disappears, and one can get around the 
first one by working with smart algorithms and small 
enough systems.) 

The first objection is a very basic one: experimentally, 
glassy systems exhibit a non-equilibrium behavior, which 
requires a dynamical description. Quite often, they ex- 
hibit a special type of dynamical behavior called aging, 
i.e., the property that extensive one-time quantities like 
the energy, magnetization etc., are asymptotically close 
to time-independent values, whereas two-time quantities, 
like the autocorrelation functions and their associated 
linear response functions, continue to depend on the time 
elapsed after the quench even for long times. Aging, de- 
fined in this way, appears in mean-field spin glasses and 
has been exhibited in spin-glass experiments We 
will not discuss systems undergoing "stabilization" j^] 
(sometimes called "physical aging" 0), where one-time 
quantities cannot be considered close to their asymptotic 
values during typical experiments. (This phenomenon 
has been recently observed in a lattice gas model with 
constrained dynamics Q.) In aging dynamics the usual 
equilibrium properties do not hold. The analysis of some 
spin-glass mean-field models |J has suggested in partic- 
ular that the usual fluctuation-dissipation relation be- 
tween the correlation and the response should be modi- 
fied in a well-defined way. This modification, which holds 
when both the age of the system and the measurement 
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time are large, involves the rescaling of the temperature 
by a "fluctuation dissipation ratio" (fdr) , which depends 
on the relation between the two times involved E^ , , 
This fdr can be found experimentally by simultaneous 
measurements of the noise and the response on various 
time scales and age scales. 

The aim of this paper is twofold. We shall first 
show that, in finite dimensional systems with short 
range interactions, there exists an identity relating the — 
experimentally accessible — FDR to an equilibrium order 
parameter function. This static order parameter function 
is an interesting new object. We shall then discuss its re- 
lationship to the usual distribution of overlaps. Our argu- 
ment relies on a perturbation of the original hamiltonian 
by the addition of some weak — but thermodynamic — 
random perturbations. This method has been recently 
used to derive interesting properties of the overlap dis- 
tribution at equilibrium 

We use the language of magnetic systems, and denote 
by S x the spin at a point a; of a lattice of size L d in d 
dimensions. We work with classical spins which are real 
variables in a double well potential, and the Ising limit 
will often be considered for simplicity. We call H(S) 
the Hamiltonian. Our argument is rather general and we 
do not have to specify much the Hamiltonian: it contains 
short range interactions, in a d dimensional space; it may 
contain quenched disorder or not. The evolution of the 
spin dynamics is governed by the Langevin equation 
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where r] x is a white noise of variance (r)x(t)r]y(t')) — 
2TS xy S(t — t'). (We denote by angular brackets thermal 
averages, i.e., either, in the dynamic framework, the av- 
erage with respect to the realization of the random noise, 
or, in the static framework, the average with respect to 
the Gibbs measure). The system starts at time t = from 
a random initial condition. Important quantities are the 
correlation function, C{t,t') = {l/N)J2 x (S x (t)S x {t')), 
and the response function, which measures the response 
of the spins at time t to an instantaneous field at time t': 



R(t,t') = ^ 



S(S x (t)) 



5n x {t>) 
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The quantity which is measured experimentally (ther- 
moremanent magnetization) is the integrated response 

function, defined by x(M') = T J*' dt"R(t,t"). 

The fdr X(q) is obtained by considering the infinite- 
time limit of the response function, fixing the correlation 
function C(t,t') to a given value q |£-11|: 

d X (t,t') ldC(t,t') 



X(q) 



lim 

C(t,t')=? 



dt' 



df 
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The usual equilibrium dynamics is obtained by sending 
the two times t 1 1! to infinity while keeping their difference 



t = t — t' fixed. Then the correlation and response func- 
tions, C(t,t') and R(t,t') reach their equilibrium values, 
c(t) and r(r). In short-range systems this regime re- 
lates to the property of "local equilibrium", i.e., to the 
fact that any finite region of space reaches equilibrium 
locally. The Edwards- Anderson order parameter is de- 
fined dynamically by qEA = lim T ^oo c(r), and the usual 
fluctuation-dissipation theorem asserts that, for q > qEA, 
X(q) = 1 . The aging regime concerns systems with weak 
crgodicity breaking, such that the correlation C(t,t') re- 
laxes below qEA when t — > oo (at fixed t') |L5|. Then the 
FDR X(q) can become different from unity in the regime 
9 < Qea- A simple way to measure it is through a para- 
metric plot of the integrated response function versus the 
correlation ||. The ratio TjX(q) can be interpreted as 
an effective temperature |19| . 

We wish to relate the fdr to an equilibrium order pa- 
rameter. Let us add to the original Hamiltonian a per- 
turbation of the form ei?2, with 



H 2 — ^ h x S x S- 



T{x), 



(4) 



where the h x s are independent gaussian random vari- 
ables of variance one, and T is a translation of length 
L/2 in a fixed direction e, say the x axis (so that 
T(x) = x+(L/2)e). The thermal expectation value of the 
perturbation (H2), is a contribution to the internal en- 
ergy of the system which is extensive and self-averaging, 
i.e., independent (in the thermodynamical limit) of the 
particular realization of the disorder contained in either 
H or Hi. The interaction H2, which looks long-range, is 
in fact a local perturbation in a different space. Let us 
divide the space into two halves (Si and S r ) and rename 
the spins in the right hand part so that if x £ Si then 
T(x) e S r and St{ x ) = S' x . The total Hamiltonian can 
now be written as 

H(S, S') = Hi(S) + H r (S') + B(S, S') + e^2 h x S x S' x . 



xeSi 
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The hamiltonians Hi and H r refer respectively to the 
spins in Si and S r . The term B(S, S') is a surface term 
whose presence does not affect the average of H2. Drop- 
ping it, the Hamiltonian (^|) characterizes a spin system 
of size L d /2, with two spins S X ,S' X on each site, and a 
purely local interaction. Notice that in the case of disor- 
dered systems the spin systems S and S' taken individu- 
ally contain two independent realizations of the disorder. 

Since the perturbation H2 is a sum of local terms, the 
thermal expectation value (for almost all realizations of 
the disorder) (H2{t)) measured in the dynamics has a 
long-time limit which is equal to its equilibrium expecta- 
tion value. The proof of this fact is standard for systems 
with short-range interactions. We first notice that the 
free energy density f(t) must reach, at long times, its 
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equilibrium value f eq : if it were to converge to a value 
/(oo) larger than the equilibrium one, one could always 
nucleate a bubble of radius r with the equilibrium free 
energy, with a cost of free energy less than or equal to 
cr d-i _|_ (-y^ _ ^(qo))^ which becomes negative for large 
enough r. Therefore the free energy reaches equilibrium, 
as well as its derivative with respect to e, proving the 
convergence of (H 2 (t)). (Notice that we do not discuss 
here the time scale for reaching this equilibrium, which 
may become very long in some systems: what matters 
here is that it is finite when L — > oo.) 

We now compute the expectation value of the pertur- 
bation H 2 in the dynamics and in the statics. Since 
(H 2 (t)) is self-averaging, it is equal to its average over 
the random field h and all other possible quenched disor- 
der in the system, which we denote by Eh (H 2 (t)). In the 
dynamical framework, starting from the Langevin equa- 
tion in presence of the perturbation eH 2 , we express the 
average of H 2 in the Martin-Siggia-Rose formalism |2l| l 
as a path integral: 

(H 2 (t))=E h (H 2 (t)) 



= E h 



V(S)V(S) e^^YhMt^Twit), (6) 



with the dynamical action 
I[S,iS] = J dt' Y^xit') 



■ OH dH 2 - 

ob x ob x 
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Integrating by parts over the h x 's, and observing that 
the insertion of iS x (t') acts as the derivative with respect 
to an impulsive magnetic field at site x and at time t': 
5/Sh x (t'), we obtain |1 



E h {H 2 (t)) =2eJ2E h 



5h T{x) {t) 



(S x {t)S x {t')S T(x) (t')). 



(8) 



In the linear response regime j3e <C 1, the average of the 
product on far away sites factorizes up to terms of order 
e and one has 

E h x , 6 u A S x (t)S x (t')S T{x) (t')) = C(t,t')R(t,t') + 0(e) 

(9) 

Assuming that the bound holds uniformly in time (re- 
member that the large volume limit is taken before the 
large time limit) and substituting the definition (|^) of the 
fdr we obtain for large values of t 



2e(3N J dqX e {q)q = e(3N (l - J dq 
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where we have assumed lim^oo C(t, 0) = for simplicity. 
We have denoted by X e the FDR of the system with the 
perturbed Hamiltonian. Notice that this is a very general 
result that holds for every sample in the case where there 
is quenched disorder. 

We now turn to the statics. The thermal equilibrium 
average of H 2 is self-averaging with respect to disorder. 
We can thus evaluate it as follows: 



-P[H(S)+eH 2 (S)] 



/,h x S x S T ( x ). (11) 



Integrating by parts over the /i x 's, we obtain 

(H 2 ) = f3eNE h (l - 1 ^2{S x S nx) ) 2 ^j (12) 

Invoking again linear response for small e, together with 
the fact that x and T(x) are infinitely far apart in the 
thermodynamic limit, we can write 



Eh(S x S T ( x )) 2 — Eh(S x S y 



0(6) 



(13) 



where x and y are two far away spins not directly coupled 
in H 2 - We obtain then (up to higher orders in e): 



(H 2 ) = !3eNE h (l- J dq P t 



(a) i' 



(14) 



The last equality, involving the overlap distribution P e (q) 
for the perturbed system, results from the decomposition 
of the Gibbs measure into a sum of pure states charac- 
terized by a clustering propert y p[ . 

Comparing the two results, ( |l0| ) and ([l4|), for the dy- 
namics and the statics, we see that the second moments of 
the dynamical order parameter function dX e (q) / dq and 
of the static one P e (q) coincide for the system in presence 
of the perturbation eH 2 . It is straightforward to gener- 
alize this derivation to perturbations of the type H p = 
J2 x h x S x S Tl ( x ) ■ ■ ■ SVp^x), where T k (x) = x + (k/pL)e. 
(For p = 1 the perturbation is nothing but a small ran- 
dom field term.) This shows that the p-th moments of 
the two functions dX e {q)/dq and P e {q) coincide. 

Let us now consider the functions X{q) = lim^o X e (q) 
and P(q) = lim c ^o P^il)- (To be precise, we need to in- 
troduce simultaneously all the perturbations with arbi- 
trary p and strength e p , and send all the e p 's to 0.) These 
are two characteristic functions of our problem. One de- 
scribes the violation of the fdt in the out of equilibrium 
dynamics, and the other describes some equilibrium cor- 
relations. These two functions are equal, and thus an 
unexpected link between statics and dynamics is estab- 
lished. 

We now discuss the relationship between the new func- 
tions dX / dq, P(q) and the more conventional definitions 
of the fdr and the overlap distribution. Let us first 
consider the equilibrium distribution P(q). Clearly, in 
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a situation with ergodicity breaking and several nearly 
degenerate pure states, the effect of the e perturbation 
which scales as ey/L^ induces a reshuffling of the weights 
of the states. A simple example appears when there is 
an exact degeneracy due to a symmetry. For instance, 
in the case of an Ising Hamiltonian quadratic in the spin 
variables, the symmetry by reversal of all the spins im- 
plies that the overlap distribution P(q) of the unper- 
turbed system is symmetric: P(q) = P(—q), each pure 
state appearing with the same weight as its symmetric 
one in the Gibbs measure. This symmetry will be vi- 
olated by the perturbation terms H p with odd p, lead- 
ing (in the absence of further reshuffling) to the rela- 
tion lim e _ >0 + P e (q) = Po+(q) = 28(q)P(q). Suppose now 
that care has been taken of all the symmetries of the 
system, by defining a modified P(q) measuring the dis- 
tance between orbits of the symmetry group. It is rea- 
sonable to believe that, for a large class of systems, the 
reshuffling of the weights will only lift the degeneracy, so 
that P(q) = P(q). We shall refer to such systems as be- 
ing stochastically stable. Mean-field spin glasses fall into 
this category, as well as Ising ferromagnets in dimensions 
d>2 p^] , but we do not know how to characterize this 
class in general. Turning to the case of dynamics, it is 
trivial to show that the limit of e — » is smooth when it 
is taken before the limit of large times (the infinite vol- 
ume limit is always taken first). If the limits commute, 
then the static X(q) = fjj dq' P(q') is identical to the 
dynamical (fdr) X(q) of the unperturbed system, mea- 
sured with random initial spin configurations. This re- 
sults holds for Shcrrington-Kirkpatrick spin glasses, but 
is violated, e.g., in p-spin spherical spin-glass models, 
where the dynamics is dominated by infinitely long-lived 
metastable states. 

Summarizing, we have introduced a new order param- 
eter function for systems at equilibrium, which can be 
related in general, in finite dimensional systems, to the 
FDR of a weakly perturbed system. This new order pa- 
rameter is interesting since its moments are obtained as 
expectation values of extensive quantitites. It is thus 
much more robust than the usual overlap distribution, 
and will not have the same chaotic behaviour under weak 
perturbations. It is interesting to notice that in mean 
field spin glasses, the order parameter appearing from 
the replica computation is naturally related to this new 
order parameter. For stochastically stable systems, this 
new order parameter is equal to the overlap distribution 
of symmetry classes of the states, and can be measured 
experimentally. The relation implies that in any finite- 
dimensional system replica-symmetry breaking and aging 
in the response functions either appear together, or do 
not appear at all. Although we have used the language 
of the magnetic systems, the arguments put forward here 
can be generalized to systems of different nature, using 
different random perturbations of the Hamiltonian. 
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